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ABSTRACT 

We study Modified Gravity (MG) theories by modelling the redshifted matter power spectrum in a 
spherical Fourier-Bessel (sFB) basis. We use a fully non-linear description of the real-space mat¬ 
ter power-spectrum and include the lowest-order redshift-space correction (Kaiser effect), taking 
into account some additional non-linear contributions. Ignoring relativistic corrections, which are 
not expected to play an important role for a shallow survey, we analyse two different modified 
gravity scenarios, namely the generalised Dilaton scalar-tensor theories and the f{R) models in 
the large curvature regime. We compute the 3D power spectrum C|(/ci, ^ 2 ) for various such MG 
theories with and without redshift space distortions, assuming precise knowledge of background 
cosmological parameters. Using an all-sky spectroscopic survey with Gaussian selection function 
ip{r) oc exp(—r^/rg), ro = 150/i“^Mpc, and number density of galaxies N = 10“"^ Mpc“^, we 
use a analysis, and find that the lower-order < 25) multipoles of Cf{k, k') (with radial modes 
restricted to fc < 0.2/iMpc“^) can constraint the parameter at a level of 2 x 10“^(3 x 10“®) with 
3(7 confidence for n = 1(2). Combining constraints from higher £ > 25 modes can further reduce 
the error bars and thus in principle make cosmological gravity constraints competitive with solar sys¬ 
tem tests. However this will require an accurate modelling of non-linear redshift space distortions. 
Using a tomographic P{a)-m{a) parameterization we also derive constraints on specific parameters 
describing the Dilaton models of modified gravity. 

Key words: ; Cosmology-Modified Gravity Theories - Methods: analytical, statistical, numerical 


1 INTRODUCTION 


The apparent accelerated expansion of the Universe jPerlmutter et al.iri999l : iRies et ailll998h can be explained within General Relativity (GR) by 
introducing a finely tuned cosmological constant. However, there are alternative e xplanations for this phenomenon, in cluding modified gravity 
theories. While the laws of gravity are not well constrained on cos mological scales l ljovce et al.ll2014l:IClifton et alj|201 3h. modification of GR ar e 


tightly constrained in the solar system JOstriker & Steinhardtll2003) or at astrophysical scales i 

Jain. Vikram & Sakesteinll2013l:IVikram et al.ll2013h. 

As shown by various authors (lBertschingeill2006l;ISong. Hu & Sawickill2007l;lBrax et al. 

I 2 OO 8 I) the background dynamics in various dark en 


ergy and modified gravity models are nearly indistinguishable. Thus it is important to investigate the evolution of perturbations in these mod¬ 
els. The studies of perturbation theory in modified gravity models, in principle, can be classified in two diffe r ent frameworks: the parametri c 
approach and the non-parametric method, e.g. the principal component analysis jZhao et al.ll2008L 120091 l^ld : lHoii^l201 ll : iHall et alj|2013h . 
Several parametrizations of modified gravity have been proposed for the evolution of linear perturbations. Such parametrizations typically in¬ 
volve two functions v(k, a) and 7 (fc, a) that both depend on the scale factor a and wave number k, when they describe the modification to 
the Poisson equation for the metric potentials. We will use instead the recently proposed parametrization of the coupling to matter /3(a) and 
the mass m(a) of the scalar field, which only depend on the scale factor (but give rise to both time and scale dependences in the resulting 
modification to the Poisson equation). The 0(a) — m(a) offers a unifie d approach to study various modified gravity models including f(R) 
gravity and the symmetron and Dilaton models l lBrax. Davis & Lill2012h . Parametrization based on effective field theory (EFT) has also been 
considered jCheung et alj 1200^1 . Besides the recent progress on the construct i on of parametrizations, m any observational windows have re¬ 
cently been proposed, such as the galaxy clustering jPogosian & Silvestrill2008| : lOvazu. Lima & HiJbOOSh . Integrated Sachs-Wolfe (ISW) effect 
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in Cos mic Microwave Background (CMB) anisot ropies ( IZhang|l2006h. the ^alaxv-I SW cross correlation dSong, Pe iris & Hu| 2007|). cluster abun 


dance jjain & Zhand 200^ Loinbriseret_^ 2012 ), peculiar velocity Li_et_^ 2000l). redshift-space distortions jjennings et_al 


_ __ _ _ __ _ _ _ _ _ ^_^j[| 201 ^ Gi^ zo et alj 

l2008h . weak -lensing Heavens. Kitching & \ferdel 2007l : ISchimdtll2008 i lTsuiikawa & Tatekawall2008ti . 21cm observations (iHall et al.ll201.'lli . matter 
bispe ctrum l GiKMarirLeLalJIlOllI), et c. In addition, recently some N-body simulation algorithms in modified gravity mod els have been devel¬ 
oped fehaolEoI^ Li. Mota & Barrowll201lh . As shown by various authors jSong. Peiris & Hij|2007l : Ikombriser et al.lf2012h . at the WMAP reso¬ 
lution the modification effects on the CMB mainly come from the ISW effect, which becomes prominent on the super-horizon scales. However, 
due to the unavoidable cosmic variance on large scales, the constraints from these effects are not significant. On the other hand, since the typical 
modification scales are in the sub-horizon range, several studies show that the most stringent cosmological c onstraints come from the large -scale 
structure data sets. For instance, in the case of f{R) gravity one obtains |/ho| < 6.5 x 10“® , 95%C.L. JPosset. Hu & ParinsonI 120140 from 
the combined analysis of the CM B temperature p ower spectrum, the galaxy power spectrum and the baryon acoustic oscillations measurements, or 
l/flol < 4.6 X 10"® ,95%C.L. jBel et al.ll20l4 from the analysis of the galaxy power spectrum through the clusterin g ratio. This is competitive 
with solar system constra ints but astrophysical constraints in dwarf galaxies provide tighter bounds | < 5 x 10“^ djain. Vikram & SakesteinI 
l2013l:IVikram et alj2013h 

The past few decades have seen a rapid progress in large scale galaxy surveys. SDS^and 2dFGS!flopened a new horizon in modem cosmology 
by mapping three-dimensional positions of millions of galaxy. The BOSS and DB S represent current state-of-the art galaxy surveys, together with the 
recently completed Wiggle^ Future surveys such as the EuclicFll lLaureiisll201 ih . LSST0 and Wide-Field InfraRed Survey Telescopy will measure 
galaxy clustering with greatly increased statistical power. They will test the theory of General Relativity (GR) on cosmological scales. One way to 
do so is to examine the growth of structure using th e 7 -paramete r (to be introduced later). Previous calculations suggest that Euclid can constrain 
7 to 0.01 (for ACDM 7 corr esponds to 7 ~ 0.55) 1 Lindej 20051; Heavens. Kitching & Verde| [200% . Generic modification of gravity may require 
more than just one parameter dAmendola. Kunz & Sap on ^ 200 djFemeira^^kotdi^ 201 ^ but previous studies have shown that a Euclid-type survey 
could measure these parameters to high precision dPaniel et al. 2010l : Amendola et alj201 lb . 

The power spectrum of density inhomogeneities in the nearby Universe and the temperature fluctuations in the cosmic microwave background 
(CMB) sky carry the bulk of the cosmological information and are routinely employed in analysing data from CMB experiments as well as large- 
scale structure (LSS) surveys.They carry complementary cosmological information; all-sky CMB observations such as NASA’s WMAF(3|and ESA’s 
Plancl(f| experiments primarily probe the distribution of matter and radiation at redshift z = 1300. Such surveys are projected or 2D surveys. The 
large-scale surveys such as those obtained with ESA’s Euclid mission will give us a window at a lower redshift range z ~ 0 — 2 and are designed 
to provide a three-dimensional (3D) map of the Universe using more than fift y million galaxies w ith spectroscopic redshifts. Ther e are other LSS 
surve ys that are either being planned or in various stage s of developments e .g. jRassat et alJiTOOSh for Dark Energy Survey (DES j^Tvson & LSSlj 
l2004h for Large Synoptic Survey Telescope (LSST) and dSchlegel et al.l200% for Baryon Oscillation Spectroscopic Survey (BOSSl’^^i These surveys 
are designed to map the local dark Universe using weak-lensing of galaxies as well as studying baryonic oscillation features in the matter power 
spectrum. This has motivated in recent years a flurry of activity in dev e loping 3D power spectrum analysis f or cosmological data using spherical 
Fourie r -Bessel (sFB) decomposition ( see e. g. Heavens_&_Ta^lm (1995); 

20K); [shapiro. Crittenden. Percivall ( |20 1 ih ; iRassat & Refregiei 2012 ); 

20 131) ). The sFB power spectrum recovers the 3D power spectrum at each wave-number k and its angular dependence is encoded in the angular 
momentum £. _ 

In this paper we use the results derived in dPratten & Munshill2013h . to derive the 3D density power spectrum for models of modified gravity 
theories. The primary goal is to check to what extent such power spectrum analysis can be used to constrain departure from GR using future redshift 
surveys. 

This paper is organised as follows: in ij^we review different types of modified gravity theories with special attention to f{R) gravity theories 
as well as the models known as the generalised Dilaton models. We introduce the 3D power-spectrum in s sFB transform, which provides 

a natural framework for treating redshift space distortions. For the specific baseline ACDM cosmology we take the following parameter values: 
h = 0.73, Dm = 0.24, Dde = 0.76, Dk = 0, iude = —1, rrg = 0.76, ris = 0.958. 


Castro. Heavens & Kitching 1 2005h: Erdogue et al. ll200d); 


Lanuss. Rassat & Strack 1 2012h : lA 


Abramo et al. 


srorev et al 


ilJ (120121)!^ Pratten & Munshi 


^ http://www.sdss.org/ 

^ http://www.roe.ac.uk/ jap/2dty 
® http://wigglez.swin.edu.au/site/ 
http://sci.esa.int/euclid 
® http://www.lsst.org/lsst/ 

® http://wfirst.gsfc.nasa.gov/ 

^ http://map.gsfc.nasa.gov/ 

® http://www.rssd.esa.int/planck 
® http://www.darkenergysurvey.org/ 
http://cosmology.lbl.gov/BOSS/ 
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Probing Modified Gravity Theories 3 


Table 1. Paiameters describing the Dilaton Models considered in our study. The parameters are used to define the scalar potential V (</?) and the coupling function 
A{ip) through the {/3(a), m(a)} parameterization defined in Eqs.^5)-(8). 


Model 

mo [/iMpc 

r 

/3o 

s 

(A1,A2,A3) 

(0.334, 0.334, 0.334) 

(1.0,1.0,1.0) 

(0.5, 0.5, 0.5) 

(0.6,0.24,0.12) 

(B1,B3,B4) 

(0.334,0.334,0.334) 

0-0,1.0,1.0) 

(0.25,0.75,1.0) 

(0.24,0.24,0.24) 

(C1,C3,C4) 

(0.334,0.334,0.334) 

(1.33,0.67,0.4) 

(0.5, 0.5, 0.5) 

(0.24,0.24,0.24) 

(E1,E3,E4) 

(0.667,0.167,0.111) 

(1.0,1.0,1.0) 

(0.5, 0.5, 0.5) 

(0.24,0.24,0.24) 


2 MODIFIED GRAVITY THEORIES 

Henceforth, we consider modified-gravity theories that differ from the standard A-CDM cosmology by adding a new degree of freedom, which 
can be associated to a scalar field ip. This new component can drive the accelerated expansion of the Universe at late times, typically through the 
nonzero value of the minimum of its potential, although this can be seen as introducing a simple cosmological constant. The main difference from 
the A-CDM cosmology and quintessence scenarios, which only modify the background expansion, is that the fifth force induced by the scalar field 
also modifies the growth of structures, typically giving a new scale dependence to the linear growing modes of density perturbations and accelerating 
the collapse of overdense regions. 

These modified-gravity theories can be classified in three broad categories: Chameleon, K-mouflage and Vainshtein scenarios, according to 
their nonlinear screening mechanism that ensures convergence to General Relativity in small-scale and high-density environments such as the Solar 
System (as the high-accuracy Solar-System measurements provide tight constraints on local deviations from General Relativity). The K-mouflage 
and Vainshtein mechanisms rely on non-standard kinetic terms that drive spatial gradients of the scalar field to zero in high-density regions and 
suppress the fifth force. Chameleon scenarios typically contain additional couplings between the scalar field and the metric, or new geometric terms 
beyond the Einstein-Hilbert gravit ational action. These two equivalent des criptions can be captured by the {/3(a), m(a)} parameterisation which 
will be used throughout this work ( iBrax, Davis & Lil2012l:lBraxetal.ll2012h . 

In this paper, we focus on two different models of the Chameleon type. In the Dilaton models, the coupling between the scalar field and the 
metric depends on the scalar field value, so that in dense environ ments where the scalar field is driven to zero, the coupling also vanishes, which 
suppresses the fifth force as in the Damour-Polyakov mechanism jPamour & Polvakovlll994) . In the f(R) models, the coupling of the scalar field 
to matter is constant but its effective potential depends on the environment and its mass becomes large in high-density regions. This suppresses the 
magnitude of the fifth force through a Yukawa screening. In all chameleon cases, these modifications of gravity induce a global enhancement of 
the effective force of gravity, due to the fifth force, which directly translates into an increase of structure formation. In this Section, we review the 
Dilaton and the f{R) models and we describe parameterization in the context of large-scale structure formation. 


2.1 Gravity in Dilaton models 

The Dilaton theories of modified gravity are chameleon models with the the Damour-Polyakov property dPamour & Polyakov! 1 1994 h. according 
to which the coupling between the scalar field ip and the rest of the matter components approaches zero in dense environments ( |Pietronill2005l : 
lOlive & PosDelovl2008l : lHinterbichler & KhourvI 2010h . In contrast with the f{R) theories described in Sec. l2.2l below. the scalar field here takes on 
a small mass everywhere and thus mediates a long-range (screened) force. These Dilaton models are scalar-tensor theories, where the action defining 
the system takes the general form 


^ ~ J d*xy/—g 


^R-l{Vpf-V{p)-At 


-b 




( 1 ) 


where Mpi = (SttC/n) is the reduced Planck mass (in natural units), g is the determinant of the Einstein-frame metric tensor and g the 
determinant of the Jordan-frame metric tensor which is given by the conformal rescaling 

gfiv = A {p>)gfi,v ( 2 ) 


The various matter fields tpm are governed by the Jordan-frame Lagrangian density and the scalar field ip by the Einstein-frame Lagrangian 
density = —1/2(V</?)^ — V(p), with the scalar-field potential V{p). There is no explicit coupling between matter and the scalar field and the 
fifth force on matter particles due to the scalar field arises from the conformal transformation ^ (more precisely, through gradients of A). In the 
Lagrangian presented in Eq. (D we explicitly wrote the cosmological constant term Aq, so that the minimum of V{p) is zero, which is reached for 
p —>■ oo, but this could also be interpreted as the non-zero minimum of the scalar field potential. 
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In the original Dilaton model, the potential of the scalar field V((p) and its coupling A(ip) with the metric have the following functional forms: 


V(p) 


Vo exp 



A((p) = 1 + 


A2 pP' 


(3) 

(4) 


where {Vo,A 2 } are the two free parameters. In dense regions where « 0, the coupling to matter is negligible, and gravity converges to GR. 
However, the field nevertheless mediates a long range gravitational force that has an effect elsewhere, in less dense environments. This model 
can be generalized to a greater class of Dilaton models, by keeping the coupling function as in Eq. 0 but considering more general potentials. 
Then, instead of specify i ng the model by its pote ntial V{(p) it is convenient to define the model by the tomographic parametrization {/3(a), m{a)} 
terax, Davis & Liir2012l : ISrax & Valageasll2013h . in terms of the scale factor a(t), where the coupling /3(a) and the scalar field mass m(a) are 
defined as 


/3(a) = /3[^(a)] = Mpi^(^), 

(5) 

m^(a) = m^lip(a),p(a)] =+ ■ 

(6) 


In this paper we consider the simple forms 


m(a) = mo a 


/3(a) = /3o exp 


- 1 


3- 2r 


(7) 


with 

9A2flmoHo ,o^ 

S o 9 * 

[The exponential potential Eq.0 corresponds to r = 3/2.] The values of the f ree parameters | mp , r, / lo , a | that enter Eg. (fTt are displayed in Table 
[3 The models {A,B,C} were chosen such as to correspond to those studied in iBrax & ValageasI ( 1201 3tl and iBrax et al.l ( l2012h . where detailed com¬ 
parisons between numerical and analytical calculations are presented. T he models |A,B,C,E| pr obe the dependence on {s, /3o, r, mo} respectively, 
other parameters being fixed (instead of the models “D” considered in iBrax & Valage^ ( 12013h . which probe the dependence on mo at fixed A 2 , 
we introduced the models E that probe the dependence on the parameter mo at fixed {s, /3o, r}). These models probe deviations from the LCDM 
cosmology of less than 20%, in terms of the matter power spectrum. 

In these Dilaton models, the coupling function A is always very close to unity, so that most Einstein-frame and Jordan-frame quantities (e.g., 
Hubble expansion rates or densities) are almost identical. Indeed, from Eqs. 0, 0 and 0 we obtain A ~ 1 -|- /3^/(2A2) and A 2 ~ (cmo/iTo)^- 
Solar System tests of gravity imply that mo > IO^Hq/c, whence A 2 > 10®, H/cm < 10“® and |A — 1| < 10“®. Therefore, the Jordan-frame and 
Einstein-frame scale factors and background matter densities, related by a = Aa and p = A~'^p, can be considered equal, as well as the cosmic 
times and Hubble expansion rates. (However, in this section we work in the Einstein frame, where the analysis of the gravitational dynamics is 
simpler.) 

In the Einstein frame, the Eriedmann equation takes the usual form, 

3MpiH = p ptp Pa, (9) 


where we consider the matter and scalar field components and the cosmological constant contribution pA- One can check that the scalar field energy 
density is negligible as compared with the matter density, p,plp ^ 10~®, so that the Friedmann equation 0 is governed by the matter density and 
the cosmological constant and we recover the A-CDM cosmological expansion, 3MpiH^ = P + pA, up to an accuracy of 10~®. In the Newtonian 
gauge, the perturbed metric can be written as 

ds^ = —(1 -I- 2^)dP + a^(t)(l — 2^)5ijdx'dx\ (10) 

where $ and are the Einstein-frame metric gravitational potentials. One can check that the impact of the scalar field fluctuations on the metric 
potentials are again negligible, as |(5p^9|/|5p| < 10“®, and we have within a 10“® accuracy $ = where q^N is the Newtonian potential 

given by the Poisson equation, 

-^q^N = 47rC/NOp = —^^3-0, (11) 

where S = dp/p is the matter density contrast. However, the dynamics of matter particles is modified by the scalar field, which gives rise to a fifth 
force given by Fa = —PV In A, that is, in the Euler equation we must add to the Newtonian potential q/x a fifth-force potential q/A = In A 
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Probing Modified Gravity Theories 5 


that is not negligible and can lead to 10% deviations to the matter density power spectrum for the parameters given in Table [T] (indeed, whereas 
IA — 1| < 10“® is negligible as compared with unity, it is not negligible as compared with |T'n|/c^ < 10 “®). 


2.2 Gravity in f{R) theories 

In models known as the ffi?) gravity, the Einstein-Hilbert action Savt is modified bv promoting the Ricci scalar B. to a function of i?, JBuchdahll 

I I I II I i‘ ‘ ' * * II I ^ 


1 1 9701: IStarobinskvIl 1 980LI 200 I : IHu & Sawickill2007h . The new action S for the f{R) gravity theories can be written as: 


S = 


J <Tx^f^ 






( 12 ) 


where we explicitly wrote the cosmological constant contribution Aq, although it is often included within the function f{R) [with our choice 
f{R) describes the deviations from GR and from the ACDM cosmology]. In this section, contrary to the previous section I tTI where we studied 
Dilaton models, we denote with a tilde E instein-frame quantities instead of Jordan-frame ones, because we now work in the Jordan frame. In the 
parameterization of lHu & Sawickil ( l2007h . the functional form f{R) can be expressed in the high curvature limit as 


f{R) 


fHo f ^ 

n R^ ’ dR J^ORn+i- 


(13) 


The two independent parameters, fiig < 0 and n > 0, can be constrained by observati ons. In the above ex pression, Ro is the present value of the 
Ricci scalar for the cosmological background. Note that this parametrization and that of IStarobinskvl l l2007h both reproduce the same results in the 
large curvature regime. 

The f{R) theories of gravity invoke the Chameleon mechanism to screen modifications of GR in dense environments such as in our Solar 
System. In this model, this occurs by requiring that all extra terms vanish in high curvature environment, i.e by requiring f{\R\ ':$> | J?o|) —>■ 0. In 
such a theory, the background expansion follows the A-CDM dynamics (as observational constraints imply \fRg \ <C 1), and the growth of structure 
is only affected on intermediate and quasi-linear scales. 

There is an essential connection between the formulation of the f{R) theory presented above, and scalar-tensor theories of modified gravity. 
Upon the coordinate rescaling = A~^ (recall that in this section is the Einstein-frame metric), with A(}p) = exp[/3v5/Mpi] and 

(3 = 1/x/G, the f{R) modifications to GR are re-cast as arising from contributions of an extra scalar field i/p, that is subject to a potential V (}p) given 
by: 


V{V>) 


Ml, ( RfR-f{R) \ 

2 V G + fR? )’ 


fn = exp 


Mpi 


- 1 . 


(14) 


In that sense, f{R) theories are equivalent to a scalar-tensor theory expressed in the Einstein frame dChiballTOO.'ll : iNunez & Solgnaikll2004h . In this 
new formulation, the screening mechanism takes another form: the mass of the scalar field grows with matter density, and a Yukawa-like potential 
suppresses the fifth force in dense environments. This can be conveniently reformulated by saying that screening takes place wherever the scalar 
field is small compared to the ambient Newtonian potential. 

It turns out that all chameleon-like models such as f{R) theories can again be parameterised by the value of the mass m{a) and the coupling 
/3(a) of the scalar field, in terms of the scale factor a and the associated background matter density p{a) = SClmoHlMh/a^. With the specific 
functional form of f{R) given by Eg.lllSt, we can directly relate {n, fptg } to {/3(a), m(a)} via: 


m(a) = mo 


/ TSJao + flmoct 
\ TUao + 


(™+ 2)/2 


mo = 


Ho / flmO + TDaO 


c y {n + 1)|/ho| ’ 


/3(a) 


1 

7!' 


(15) 


In this paper, we consider values of n = {1,2} and \fRo\ = {10“^, 10“®, 10“®}. The larger value of |/ho| is cun'ently ruled out by other 
independent probes, so this serves as a consistency test. 

As for the Dilaton models described in Sec. 12.II the f{R) models that we consider in this paper follow very closely the A-CDM cosmology at 
the background level, because Ifug | <IC 1. Indeed, from the action dl2t one obtains the Friedmann equation as 

3M|i - Jr{H^ + H) + f/Q + JrrHR^ =P + pk, (16) 

where the dot denotes the derivative with respect to cosmic time t and fnn = f /dR^. In the background we have R = 12H^ -|- 6H and we 
can check that all extra terms in the brackets in Eg. ( 116b are of order \fRo\H^, so that we recover the A-CDM expansion, = P + PA, up 

to an accuracy of 10“"^ for \fRg\ < 10“^. Moreover, the conformal factor A{ip) is given by A = (1 -I- fR)~^^^, so that |A — 1| < 10“^ and 
the background quantities associated with the Einstein and Jordan frames can be considered equal, and equal to the A-CDM reference, up to an 
accuracy of 10““^. Considering the metric and density perturbations, we can again write the Newtonian gauge metric as in Eg. dlOb (but this is now 
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the Jordan-frame metric). Then, in the small-scale sub-horizon limit, the modified Einstein equations lead to 
2 2 

where 5fa = fn — fp and is the Newtonian gravitational potential defined as in GR by Eg. dl lb . Thus, because we work in the Jordan frame, in 
contrast with the Dilaton case presented in Sec. 1 12.11 the modification of gravity directly appears through the metric potentials Ea.(ll7t. Finally, the 
dynamics of the matter particles is given by the geodesic equation, or the Euler equation in the large-scale single-stream limit, where the Newtonian 
potential that appears in GR is replaced by the potential given in Eq. G3- 


2.3 Modified gravity and structure formation 

2.3.1 Impact of modified gravity on 3D matter clustering 

As we have seen in the previous sections, Dilatons and f{R) theories reproduce the smooth background expansion history of the standard A-CDM 
cosmology (up to an accuracy of 10“"^ or better that is sufficient for our purposes). To distinguish between, and hopefully test such competeting 
gravitational theories it is thus necessary to analyse the evolution of matter or metric perturbations. To lowest order in cosmological perturbations, 
non-standard gravitational scenarios effectively result in a time- and scale-dependent modification of the Newtonian constant (Jn- In the nonlinear 
regime, the modifications become more complex as they become sensitive to the screening mechanism, which depends nonlinearly on the envi¬ 
ronment and modifies the form of the equations of motion (e.g., the effective Poisson equation is no longer linear). These modifications induce a 
distortion of the dynamical as well as statistical properties of the matter clustering statistics. To this end, we focus on the signature of these effects on 
the matter density power spectrum P{k, z), or more precisely its expansion Ce{k) on spherical multipoles, which is well suited to the cosmological 
analysis of wide galaxy surveys. (We also include redshift-space distortions through a simple approximation that is exact at linear order and includes 
some non-linear contributions.) 

Thus, we need a way of modelling of t he matter power spectrum that applies to the standard A-CDM scenario as well as these modified-gravity 
models. We use the appr oach developed in iValageas. Nishimichi & Taruval 20131) . for the A-CDM cosmology, and next extended to these modified- 
gravity cosmologies in jBrax & Valageasll2013l) . This method combines the results from one-loop perturbation theory with that from halo model 
predictions. By construction, th is power spectrum agree s with perturbation theory up to order P|. In the A-CDM cosmology, this corresponds to 
the standard one-loop diagrams iBemardeau et alj ( |2002|) . In the modified gravity models, the linear propagators and the vertices are modified, with 
new scale and time dependences. In addition, the effective Poisson equation, obtained after integrating over the scalar field, becomes nonlinear. At 
one-loop order this gives rise to a new diagram for the power spectrum. Beyond the perturbative regime, predictions from a suitably altered halo 
model are incorporated in the high-fc limit. The impact of the modified gravity in the nonlinear dynamics is take n into account through the halo m ass 
function but the impact of the modifie d gravity on the halo profiles is ignored, i.e., we keep the NEW profile from lNavarro. Frenk & Whita ( ll996h and 
the mass-concentration relation from IValageas. Nishimichi & Taruval l l2013ll . The resulting matter density power spectrum Pssj k) has been tested 
again st numerical simulations and found to be in agreement for the entire available k range from simulations k < 3/i“^Mpc dBrax & Valageaj 
l2013h . 


2.3.2 Numerical results for the linear and non-linear 3D power spectra 

In the modified gravity scenarios, the linear growing mode D+ {k,t) and the linear growth rate f{k,t) = 9In D+ /din a become scale dependent. 
We show the linear growth rates f{k, z) in Figure-[T]for the f{R), (n = 1), models and in Figure-l^for the Dilaton models. For the Dilaton models, 
we only show those models that exhibit maximal departure from GR predictions. The modified gravity models that we consider in this paper amplify 
and accelerate the growth of large-scale structures, so that their linear growth rates are greater than the ACDM one. The relative importance of the 
fifth force also increases with time and is mostly relevant in the dark energy era, at z < 2. In addition, we can see in Figures-Q] and -[2] the scale 
dependence generated by these modified gravity models. The deviation from the ACDM growth rates decreases along with I/rqI- For the dilaton 
models the deviation decreases on these scales as B4 > A3 > E3 > Cl. 

The non-linear matter power spectra for the modified gravity theories are displayed in Figure-[3 The deviations from the ACDM power spectrum 
peak at the weakly non-linear scales k ~ 5h,Mpc“^, due to the amplification by the non-linear dynamics. Moreover, at very large scales (much 
beyond the Compton wave length) all models converge back to GR and to the ACDM cosmology. At high k,k> lO/iMpc”^, we may underestimate 
the deviation from the ACDM power spectrum because we neglect the impact of the modified gravity on the halo profiles. However, the analysis 

i iresented in thi s paper is restricted to linear or weakly non-linear scales, k < 0.2hMpc“^, where the detailed shape of halo profiles plays no role 
Valageasll2013l) and our modelling of the matter density power spectrum is reliable. In agreement with Figure-[T] the deviations from the ACDM 
power spectrum increase with \fRg \ and in the order B4 > A3 > E3 > Cl. For the f{R) theories the deviations are smaller for a greater exponent 
n, because this gives a faster increase of the scalar field mass m at higher redshift, as seen in Eg. ( 115b . whence a faster convergence to ACDM. 
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J{k,z)/kR{k.:) |/o| = 10-* n = l 



k [/?.Mpc *] 


I{k.z)/fcR{k.z) |/o| = 10-' n = l 



k [/? Mpc *] 


J{k,z)/fGR{k.z) l/ol = 10 n = l 



k [/?.Mpc *] 


Figure 1. We plot the linear growth rate f(k^z) = d\n Dj^(k, z)/d\na, normalised to the A-CDM value, for /(R) models as a function of the wave number 
(in units of hMpc~^) and of the redshift 2 . The functional form for the f{R) model is described in Eq. (H, with n = 1. The parameter is fixed at —10 
(left-panel), —10“^ (middle panel) and —10~® (right-panel), respectively. 


f{k,z)/fGn{k,z) A3 f{k,z)/fGR{k,z) B4 



Figure 2. We plot the linear growth rate f{k, z) = 5 In D+ (fc, 2 )/^ In a, normalised to the A-CDM value, as in Figure^but for the Dilaton models. The top panels 
correspond to the Dilaton models A3 (left) and B4 (right) models. The bottom panels correspond to the Dilaton models Cl (left) and E3 (right) respectively (i.e. only 
the extreme case from each class of models is shown). 
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Figure 3. The ratio of the matter power spectra in modified gravity theories and in the ACDM model is plotted as a function of the wave number k. The left panel 
shows the results for f{R) gravity and the right panel for Dilaton gravity models. The redshift in both panels is fixed at 2 : = 0.5. We show two different models of 
/(/?) gravity that correspond to n = 1 (solid-lines) and n = 2 (dashed-lines). For each model we consider the three different amplitudes = —10“^, —10“® and 
—10“® (from top to bottom), respectively. The parameters of the Dilaton models displayed in the right panel are given in Table-IT] Four different models are shown. 
Only the extreme version (i.e., with the greatest deviation from ACDM) of each type of model is shown. 


3 LINEAR REDSHIFT SPACE DISTORTIONS (BEYOND KAISER EFFECT) IN MODIFIED GRAVITY 

3.1 Spherical Fourier-Bessel Formalism 

3.1.1 Definition of the spherical Fourier-Bessel transform 


Spherical coordinates are often a natural choice in the analysis of cosmological data sets as they can, by an appropriate choice of coordinates, be 
used to place the observer at the origin of the analysis. As upcoming surveys promise to yield both large (i.e. wide angle) and deep (i.e. large radial 
coverage) coverage of the sky, we require a simultaneous treatment of the extended radial cover age and spherical sky geometry. A natural basis for 
such an analysis is given by the spherical Fourier -Bessel formalism. In this section, we follow jHeaven sl l2003l : ICastro, Heavens & Kitchindl2005l : 
iRassat & Refregi^l2012l : IPratten & Munshiir2013h and detail the conventions used in this paper. 

Consider a homogeneous 3D random field ip (r, fl) such that Q denotes the angular coordinate on the surface of a sphere and r denotes the 
comoving radial distance. The eigenfunctions of the Laplacian will be constructed from products of the spherical Bessel functions of the first kind 
ji{kr) and spherical harmonics with eigenvalues of —k^, for a 2-sphere. Assuming a flat background Universe, the sFB decomposition of 

the homogeneous 3D field will be given by 




POO 

/ dk Y, 'i'emik)kMkr)Ye^iCl), 

{em} 


with the inverse relation given by 


’i'emik) = \ - drrY dQ.-^{r,n)kj,{kr)Yl^{D.). 


(18) 


(19) 


This is something of a spherical polar analogue to the conventional Cartesian Fourier decomposition. In particular, defining the normalisation of the 
3D Fourier transform and power spectrum as 


'l'(r) = 


1 


dke"’' *' ^(k), (T'(k)v&*(k')} = P(fc)<5D(k - k'), 


(27r)3/2 

the sFB coefficients and the Fourier modes can be related as 


’ktm(fc) = /fc J dCl Y;^{Ci) <f{k, Cl), i 


( 20 ) 


( 21 ) 


{^7n} 


( 22 ) 


while the two power spectra obey 

C^tru{k)n'm'{k')) = Ct.{k) Soik - k') Sw Smm> 


with Ci{k) = P{k). 
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3.1.2 Finite-depth surveys 

In reality, we will often want to consider a cosmological random field that is only partially observed due to a finite survey volume. In this instance, 
we can construct the observed field 'P(r, fl), which we denote with a tilde, by multiplying the original field ^'(r, Cl) with a selection function ip{r) 

^{r,Cl) = 'i>{r,Cl)p{r). (23) 

The sFB coefficients of this finite-depth field can be related to those of the field T' by 

roo cy poo 

/ dk'Wdk,k')^t,n{k') with We{k,k') = — drr'^ifi{r)kje{kr)k'ji{k'r). 

Jo Jo 

The introduction of the selection function means that the homogeneity criterion is not valid in the radial direction and the observed sFB power 
spectrum will now given by 

poo 

{^im{k)^;,^,{k')) =Cdk,k')6ee,Srr,m', with C,(fc,fc')= / dk"Wdk,k")Wdk',k”)P{k''). (25) 

Jo 

Typically we will often be interested in the diagonal modes for which k = k' as the sFB power spectrum falls off rapidly away from the diagonal. 


(24) 


3.2 Spherical Fourier-Bessel Formalism: Applications 

As previously mentioned, the 3D approach will be particularly important for the analysis of cosmological data from future wide-field surveys. 
This is especially true for large angular scales in which the plane parallel approximation, the distant observer approximation or, equivalently, the 
high I approximat ions to the spherical harmonics are all inadequate. Conventionally, redshift space distortions (RSD) are typically studied in 3D 
jFisher et alll 19951) whereas weak-lensing has traditionally been studied with projected surveys (2D) due to a lack of, or uncertain, photometric 
redshift information about individual sources. Due to the information available, many of these studies were also limited to small patches of the sky 
and therefore invoked the flat-sky approximation. With the advent of surveys that can provide accurate photometric redshift information across wide 
areas of the sky, there has been a growing demand for techniques that can make use of this full-sky 3D information. For example, an early approach 
to incorporating this 3D information into the 2D projected surveys was to invoke some form of tomographic reconstruct ion in wh i ch the sources are 
divided up into slices at different redshifts and a 2D analysis is performed in each of these redshift bins. More recently. Illeavensl ( l2003h proposed a 
genuine 3D formalism for weak-lensing surveys based on the spherical Fourier-Bessel expans ions. These studies were later extended to a detailed 
description of weak lensing observables on the full 3D skv ICastro. Heaven s & Kitchind (l20 05h. In addition, power-spectrum estimation techniques 
have been generalised to the analysis of higher-order statistics in 3D i Munshi. Heavens & Colesll201 lal : lMunshi et al J201 Ibh . 

Another growing are a of research within the 3D approach is the cross-correlation of galaxy and weak lensing surveys with other cosmological 
observables. For example. IShaniro. Crittenden. PercivM (1201 1) studied the cross-correl ation of 3D galaxy surveys with the projected CMB in order 
to study the integrated Sachs-Wolfe (ISW) effect. Likewise, Pratten & Munshil (1201 4h detailed the cross-correlation of 3D weak-lensing with the 
projected thermal Sunyaev-ZeTdovich (tSZ) effect as a way to recover redshift information that is lost in the line-of-sight projection of the thermal 
pressure of free elections. 

The 3D approach can also be used to study and characterise bary on acoustic oscillation (BAO) features in the matter power spectrum. The sFB 
approach was first used i n this context bv |Rassat & RefregieJ ( l2012n and later extended to include linear RSD as well as non-linear effects in the 
matter power spectrum bv IPratten & Munshil 1 2013 ). 


3.3 Redshift Space Distortions 

3.3.1 Expansion to first order over peculiar velocities 


The existence of inhomogeneous structure in the Universe induces peculiar velocities that lead to distortions in the observed clustering of galaxies as 
measured in redshift space. The anisotropies generated by these distortions are known as redshift space distortions (RSD) which, together with bias 
and non-linear evolution, i nduce departu res in the measured matter power spectrum away from the configuration-space power spectrum predicted by 
linear perturbation theory ( lKaiseilll987l) . These distortions necessarily complicate the cosmological interpretation of spectroscopic galaxy surveys 
but the RSD are also one of the most promising probes for the measurement of the growth rate of structure formation and hence a useful probe for 
models of dark energy and modified theories of gravity. 

The effect of the RSD on the matter power spectrum and clustering statistics can be broadly split into two effects: the linear Kaiser effect and 
the finger of God (FoG) effect. The linear Kaiser effect is a coherent distortion of the peculiar velocity along our li ne of sight w ith an amplitude 
controlled by the growth rate. The Kaiser effect leads to an enhancement of the power spectrum amplitude at small k dKaise J 19871) . The FoG effect 
arises due to the random distribution of peculiar velocities for galaxies within virialized structure s. These peculi ar velocities lead to an incoherent 
contribution in which we have dephasing and a suppression of the clustering amplitude at high k djacksonlll972h . The effect of a peculiar velocity. 
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or departure from the Hubble flow, v(r) at r is to distort the observed comoving position in redshift space s from its true comoving position in real 
space r: 


3(r) = r + 


v(r) ■ Q. 

aH ■ 


(26) 


In the following, we denote by a superscript “s” fields that are defined in redshift space s, to distinguish them from the real-space fields. The 
redshift-space sFB transform can still be defined as in the real space case, 




ir-7 




(27) 


The conservation of matter implies [1 + 5'’(s)]ds = [1 + d(r)]dr, so that in the case of the sFB transform of the density contrast we can make the 
change of integration variable from s to r, as 

/ Q roo r j r\ noo r 

Stm{k) = — j drr^ J dD,[l + S{r,{l)] k — J ds J d{lk je{ks)Yi^(Cl). (28) 

The last term only contributes for the monopole (f = m = 0). Expanding over the peculiar velocity, we obtain 


/ 2 r°° f 

5 l^{k) = ^^-j drr^ J dn[l + S{' 




kji{kr) + ^ 7>(fcr) + ... 




Using the following perturbative expansion over powers of the peculiar velocity: 

we have 


and 


Rm}/{0,0}: 3 fJ^\k) = 5 em{k) = ^ drr^ j dQ 5{r,Cl)k je{kr)Y;^{Cl) = k J dQ S{k,n)Y;^{n) 


Sfr^\k) = ^^J^"drr^Jdn[l + S{r, j',{kr)YUCl) 


(29) 


(30) 


(31) 


(32) 


The continuity equation reads as adS/dt + V[(l + d)v] = 0. Even though (1 + 5)v is not curl-free, to obtain an order of magnitude estimate we 
may write 


which is exact at linear order, and 


.?»(!) 


(fc) - 


L ""V 


2 
TT 

,2 roo 


dn V 


,-i^ 

dt 

d5{k',n) 


{r,n)-^k^j[{kr)YUn) 


dk' k' I drr^j',{kr)j[{k'r)Yl^{h). 


(33) 


(34) 


(35) 


In the linear regime, the growth rate of the linear growing mode D+{k, t) (which usually depends on the wave number in modified-gravity scenarios) 
is defined as 


, . d\nD+ IOIxiPl , 

fik,t) = —— 

0 In a 2 0 In a ot 

In a similar fashion, we define the nonlinear growing mode Dp(k, t) and growth rate fp{k, t) from the non-linear power spectrum as 


Dp{k,t) = 




PLik,0)’ ' 2 aina ’ 

where Pp and P are the linear and non-linear power spectra, and we use in Ea.(l35t the approximation 
c)(5(k,f) 


dt 


~ fpHS, 


(36) 


(37) 


(38) 
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which is exact at linear order. This gives 


5 


s{l) 

im 


(k) 


~ -i‘k^ [ dk'k' 
TT Jo 





drr'^5{k\ Cl)fp{k', z)jt{kr)jt{k'r). 


(39) 


Thus, Eas.l l31t and J39t give the exact expression of the matter density sFB transform in redshift space at linear order. Moreover, it includes a partial 
account of nonlinear contributions. The zeroth-order term Ea.J3U over velocities includes all nonlinear contributions from gravitational clustering. 
The first-order term Ea.ll39t over velocities is exact at linear order but only approximate at nonlinear order, because of the approximations Eg. <331 
and EQ.( l38t . Nevertheless, they should capture the magnitude of nonlinear contributions to this large-scale Kaiser effect. 

Although we go beyond the usual linear-order approximation, by taking into account this nonlinear contributions, we restrict ourselves to 
large weakly nonlinear scales, because we neglected small-scale virial motions associated with the fingers-of-god effect. Thus, the redshift-space 
distortions considered in this paper correspond to the Kaiser effect associated with large-scale coherent flows. 


3.3.2 Finite-depth surveys 

In the case of finite-depth surveys, as in 1 )3.1.21 we need to multiply the density field by the selection function ip{r). We still define the selection 
function in real space rather than redshift space because it is not necessarily affected by peculiar velocities in the same fashion as the radial coordinate. 
Typically, the selection function depends on the flux and observed angular size of the objects and writing it as a function of redshift is a convenient 
approximation. Besides, it typically varies on cosmological scales, of order c/H, which are much greater than the weakly nonlinear scales that we 
aim to probe (or order 2n jk). Then, one can write ip(f) — ‘Pif, which also corresponds to neglecting logarithmic radial gradients of the selection 
function as compared with logarithmic gradients of the density field. 

Then, Eg.ll^ becomes 


{£,m}/{0,0}: dr r'^ J dCl[l + S{r,Cl)] p{r) 


kje{kr) + -|- ... 

an 


YLm 


We can again expand over the peculiar velocity as in Eo.l l30b . The zeroth-order component is given by the real-space expression l l24b . 


dk'Wf°\k,k')Spi^{k',0) with 
Jo 

9 /"oo 

Wf^\k,k') — — / dr r'^ Dp{k', z)ip{r)k ji{kr) k' ji{k'r), 
Jo 


(40) 

(41) 

(42) 


where we defined Spemik', 0) the linear density contrast today, at z = 0, and Dp{k, z) is the non-linear growing mode defined in Eo.l l37b . while 
the first-order component reads as 


dk'wf^\k,k')5L£mik',0) with wf^\k,k') = -k'^ dr r'^ Dp{k', z)fp{k', z) p{r)j'i{kr) j'fk'r). 


(43) 


We again define the auto- and cross-power spectra constructed from these harmonic coefficients as 

{K^^\k)5lf2;{k')) = cf‘^’^\k,k')6u'S^m', 

where the indices a and /3 take values {0, 1}, associated with the first two orders of the expansion over peculiar velocities. From Eos. I42b and I43b 
we obtain 


(44) 


Cf°‘'^\k,k') = j dk” Wf°‘\k,k'')wf^\k' ,k")PL{k" ,0), 

and, up to first order over peculiar velocities, the redshifted power spectrum reads as 


(45) 


C!{k,k')^Cl^°°\k,k')+cf^°\k,k')+cf°^\k,k'). (46) 

The covariance matrix C|(fc, k') is symmetric and approximately block diagonal in nature. 

The kernels are exact at linear order over the matter density fluctuations and contain some non-linear contributions. Thus, the real-space 

power is exact at full non-linear order, within the model that we use to compute the non-linear matter power spectrum, but the cross-power 

is only exact up to linear order, because it involves the approximations Eo.l l33b and Eg. I38b . To check that the non-linear contributions do not 
play a significant role, and do not degrade our predictions, we compare our results with those associated with the linear power spectrum itself. Thus, 
we introduce the kernels and Wff that source the 3D linear power spectrum Pp{k', 0) and depend on the linear growth rate D+ respectively 
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through the following convolutions: 

W"r}i\k,k') = — / drr^D+{k',z)ip{r)kji{kr)k'ji{k'r)-, (47) 

Jo 

Wl[^\k,k') = ^k'^ j drr^ D+{k',z) f(k',z)ip{r)j't,{kr)j[{k'r). (48) 

These are the linear counterparts of the kernels and defined in Eas. ll42b and d43b . The corresponding sFB linear power spectrum is 

defined through the following equation; 

Cl[°‘'^\k,k') = j dk''Wl["\k,k'')Wl[^\k\k'')PL{k'',Q). (49) 


3.3.3 Galaxy bias 

In practice, we do not observe the matter density field itself, except in weak-lensing surveys, but the galaxy distribution. Writing the galaxy density 
field as a linear function of the matter density field with a scale-independent bias, Sg (r, z) = b{z)S{r, z), while the velocity field remains unbiased, 
we recover the same expressions as in the previous section, but with the kernels given by 

Wl‘^°\k,k') = h{z)W‘^°\k,k'), Wf’‘^^\k,k’) = Wf^\k,k'), (50) 

where again the first-order term is only exact up to linear order over the density and velocity fluctuations. The same relations hold between 

the linear kernels and . Then, the expansion d46b is also a first order expansion over the ratio jJp = fp/h. 


4 RESULTS 

If we ignor e the effects introduced by the select ion function, i.e. set ip{r) = 1, and neglect redshift space distortions, i.e. set s = r, then we recover 
the result of ICastro. Heavens & Kitchind ( l2005h for the un-redshifted contributions Ci{k, k) = Pss{k). These expressions hold for surveys with all¬ 
sky coverage. In the presence of homogeneity and isotropy, the 3D power spectrum will be independent of the angular multipole i. The introduction 
of a sky mask breaks isotropy and introduces additional mode-mode couplings. The machinery for dealing with partial sky coverage and a sky mask 
is reviewed in the appendix of IPratten & Munshil ( l2013h . Note that in the above equations we neglect a number of additional non-linear terms. These 
include General Relativi stic corrections, velocity terms as well as lensing contributions. In addition, the flat sky limits of Eg. ( 145b can be found in 
dPratten & Munshill2013ll . 

For ACDM, j3{k, z) « I3{z) but in modified theories of gravity f3{k, z) = f{k, z)/b{z) is typically fc-dependant. For both ACDM and the 
modified theories of gravity considered in this paper, we use the full f3p{k, z) taken from the numerical calculations of the non-linear matter power 
spectrum. 

We choose a survey with selection function (p{r) oc exp(—r^/rg) and rg = 150/i“^Mpc. For the bias we use the model b{z) = + z. 

However, we find that our results are not very sensitive to b{z) for the survey configuration that we have considered. The results are presented for 
all-sky coverage. 


4.1 Impact of non-linear contributions and of redshift-space distortions 

We show the diagonal entries of the 3D sFB power spectrum C| {k, k), as defined in Eg.ll46b. in Figure-O for the ACDM cosmology. We also plot 
separately the contributions (equal to the real-space power) and -f (the redshift space contribution at first order over the peculiar 

velocity power spectrum). In each case we also display the linear power spectra defined by Eq. Hi. 

We find that the non-linear contributions are quite small below k < 0.2/iMpc~^ but cannot be neglected at higher k. This transition scale also 
corresponds to the harmonic i > 10. For higher £ the non-linear contribution becomes more pronounced at lower k. The redshift space distortions 
typically give a 20% contribution to the total power on the scales displayed in Figure-]?] 


4.2 Off-diagonal terms of the power spectrum matrix C| (fc, k') 

A few slices through the window functions k') and W^‘''^\k, k') that represent mode-mixing in the presence of the radial selection 

function are displayed in Figure-|5]for 1 = 2 and Figure-j^for I = 10. For higher I and k' the mode mixing is more suppressed, that is, the window 
function is more strongly peaked around k' ~ k. 
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ACDM 



k (h Mpc ^) 


ACDM 



k (h Mpc ^) 


ACDM 



k (h Mpc ^) 


Figure 4. The 3D power spectra C^{k) = for ACDM are shown as a function of the radial wave number k for various angular harmonics. 

From left to right the panels represent £ = 2,10 and 50, respectively. The topmost cui've in each panel coiTesponds to C|(fc), from Eq. j46t . Next, the two sets of 
curves, from top to bottom, correspond to the contributions from and respectively. In addition to the non-linear power C^{k) (solid-lines) 

defined by the window function <45 1 . we also show the linear power Cf^^{k) (lower-dashed lines) defined by the window function <49 F 




{k\k) 




k (h Mpc ^ ) 


k (h Mpc ^ ) 


Figure 5. The kernel functions introduced in Eq. go) are presented. The left panel shows W^^^^'^{k,k') and the right panel W^^^^\k,k') for i = 
2. Different curves with peak position shifting from left to right correspond to k' = 0.001 (solid lines), 0.025 (dashed lines), 0.1 (dot — 
dashed lines), 0.35 (dotted lines) hMpc“^. 


{k,k') 


{k',k) 




Figure 6. Same as Figure-[5]but for i = 10. 
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k (h Mpc ^) 


Cf {k,k') (k,k!) 



k (h Mpc ^) 


Figure 7. The slices of the covariance matrices are shown for f = 2. The left panel shows {k,k') as a function of k for fixed values of k'. The right panel corre¬ 
sponds to (fc, k') + (fc, k'). Various solid curves represent different fixed values of k', with k' = 0.005(solid lines), 0.05(dashed lines), 0.1(dot — 

dashed lines) and 0.35 h Mpc“^ (dotted lines) from left to right. The thick dashed lines in each panel represent the diagonal entries of the covariance matrices, 

q{k,k). 


{k,k') (k,k') (k,k') 


10^ 

£ = 10 


10^ 

£ = 10 ^ ~ T; 

l-» 1 

O ' 
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l-» 
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10^ 

- ' ' . 


l-» 

O 

_ f . 1 . 

l-» 

O 

: A. .i ; 


l-» 

o 


10 '^ 10 '^ 10 '^ 10 '^ 

k (h Mpc^^ ) fc (h Mpc^'^ ) 


Figure 8. Same as Figurel7]but for £ = 10 and with k' = 0.05(dashed lines), 0.1(dotted dashed lines) and 0.35(dotted lines) h Mpc ^ from left to right. 


The power spectra are diagonally dominated. However, a few off-diagonal terms of the matrix C| {k,k') are displayed in Figure-[7](for i = 2) 
and Figure-[8] (for £ = 10). The left and right panels show and -I- respectively. For higher k' the covariance matrix is more 

sharply peaked at fc = k'. The slices through the covariance matrix for £ = 10 as depicted in Figure-[^are more sharply peaked compared to the 
£ = 2 results for the same values of k'. These properties follow from the behaviour of the window functions FF/(fc, k') displayed in Figure-lSjand 
•[^ To compare the diagonal elements with the off-diagonal terms we have also shown C|(fc, k) in these plots (dashed-lines). 


4.3 Impact of the modified gravity models 

In Figure^^and Figure{T0]we plot the diagonal entries of the 3D sFB power spectrum C|(fc, k), as defined in Ea.ll46t. as a function of the wave 
number k for f{R) gravity theories. The results correspond to £ = 2,10 and 50 respectively. The range of fc-values probed is O.OlhMpc”^- 
0.2 hMpc“^. In Figure^^we show the results for n = 1 and in Figure{T0]we show the results for n = 2. For each value of n we show the three 
values I/rq \ = {10“^, 10“®, 10“®}. The base A-CDM model is also plotted (solid lines). In comparison to the A-CDM model all f{R) models 
have additional power at all k. We find that for all values of k and £ the redshift-space contribution -I- is positive. We will see that 

the inclusion of redshift information improves our ability to distinguish departures of MG theories from GR. The Dilaton models we consider are 
specified by the four parameters mo, r, fio and s, as given in Table-|T] The results for £ = 2, 10 and 50 are presented in Figures- fTTl H21 and 1131 
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f(R) n = l 



k (h Mpc ^) 


f(R) n = l 



k (h Mpc ^) 


f(R) n = l 



k (h Mpc ^) 


Figure 9. The 3D linear power spectra C£^(fc) that represent the diagonal elements of the 3D sFB power spectrum C£^(fc, k) as defined in Eq. i46t for f(R), n = 1, 
theories of gravity are plotted as a function of the radial wave number k for £ = 2 (left-panel), £ = 10 (middle-panel) and £ = 50 (right-panel), respectively. In each 
panel the dashed curves correspond to the three different parameter values = —10“"^, —10“^ and —10~® (from top to bottom). The ACDM result is shown 
with solid curves. Each panel displays two different sets of curves. The top and lower sets con'espond to sFB power spectra C£^(fc) with and without the redshift space 
distortions. 


f(R) n=2 


f(R) n=2 


f(R) n=2 




CO 

o 




k (h Mpc ^) k {h. Mpc ^) 

Figure 10. Same as Figurejojbut for n = 2. 



k (h Mpc 


respectively. In agreement with the 3D power spectra shown in Figure-l^ over the range k < 1/iMpc ^ and £ < 50 associated with linear and 
weakly non-linear scales, the relative deviations from the ACDM power grow at higher radial wave number k and angular harmonic 1. 


5 COVARIANCE AND 


The likelihood function C for arbitrary sets of parameter 0^ (that specify a given MG theory; e.g. {n, } in the case of f{R) theory), given the 

data vector C| (k) (which consists of the noisy 3D sFB power-spectra C| (fc) = C| (k, k)), is given by: 


£(e^|c|(fc)) 


1 

(27r)^pix/2|det C|i/2 



dk 


dk' SC!{k) CJ^{k,k') 5Ch{k') 


(51) 


Here SCI = C| ~ (Cl'^^ being the power spectrum in GR); A^pix is the size of the data vector, which depends on the angular resolution Imax 

and the number of radial bins used in the computation. The covariance matrix Cu' (k, k') is given by: 


CeedKk') = {Ct{k)Ct,{k')) - (C|(fc))(C|.(fc')> 


2 

21+1 


Cl {k, k ) + 


-T 

N 


Su'S 

m m' ? 


(52) 


where we used a Gaussian approximation. Here N represents the number density of galaxies which are assumed to be Poisson distributed. The block 
diagonal form is a result of assuming an all-sky coverage. However we found the C to be diagonally dominant. Partial sky coverage will introduce 
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Figure 11. The 3D linear power spectra {k), as in Figure[^but for the Dilaton models, are shown as a function of the radial wave number k for the angular harmonic 
^ = 2. In each panel the solid curves represent the ACDM result while the dashed curves represent the Dilaton models. From left to right the panels correspond to 
the Dilaton models A, B, C and E, respectively. The upper (lower) set of curves in each panel con'esponds to results with (without) redshift space distortions. For each 
chosen Dilaton Gravity model we plot the power spectra for the three choices of parameter values given in Table-^ The models A, B, C and E correspond to variation 
of the parameters s, r and mo respectively. 
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Figure 12. Same as FigurefTTIbut for £ = 10. 
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Figure 13. Same as FigurefTTIbut for I = 50. 
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Figure 14. The yp for the f{R) models are presented as a function of the parameter I/hq |. The left panel correspond to n = 1 and the right panel to n = 2. The two 
solid curves near the top of the panels correspond to results using non-linear power spectrum (45). The curve with higher for a given includes RSD contribution 
(nLinR) while the one (nLin) with lower doesn’t. The two dashed curves at the bottom of the panel are the results derived using linear power spectrum ED with 
and without RSD contribution denoted as LinR and Lin respectively. Inclusion of RSD improves the t)oth linear and non-linear regime. 
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Figure 15. The x^ for various generalised Dilaton models (see Table0 are shown. From left to right various panels correspond to models A, B, C and E. Each model 
corresponds to variation of a specific parameter among {s, t", mo} while keeping all other parameters fixed. The x^ in panels from left to right is shown as a 

function of s, r and mo- The line-style is the same as in Figure- Tin 


off-diagonal terms between different harmonics in the covariance matrix. A rough scaling of signal-to-noise (S/N) with the fraction of sky-coverage 
/sky is typically used in the literature: S/N oc We work with the x^ statistics defined as: 

^ /dfc / dk' SCl{k) CJ^]{k, k') SCtdk')- (53) 

w •' •' 

We assume a perfect knowledge of all background cosmological parameters and present results with and without RSD. In practice, we bin the radial 
wave number k in four logarithmic bins to avoid the covariance matrix being singular, thus replacing the integrals over k and k' by discrete sums, 
and we use the multipoles £= 1 to 25. The galaxy number density is fixed at N = 10“"^Mpc“^. We again assume all-sky coverage and the survey 
depth is fixed at ro = 150/i~^Mpc. We show our results for the f{R) theories in Figure-fTTl For the n = 1 model we find that for N = 10“^Mpc~^ 
the values of \fRg\ > 2 x 10“® can be ruled out with a 3a confidence. For n = 2 we find the constraint degrades to |/ijo| > 3 x 10“®. The 
f(R) models with higher exponent n are progressively less constrained as they converge increasingly fast to ACDM with redshift at 2 > 0. We 
show our results for the Dilaton models in Figure-fTSl We consider the problem of the estimation of each individual parameter {s, /3o, r, mo} while 
keeping others fixed, which corresponds to the model families A, B, C and E. For N = we find that the parameter values /3o > 0.8 

and mo < 0.15/iMpc“^ can be ruled out with 3cr confidence. No meaningful constraints on cr or r can be obtained at 3a level. Clearly, a joint 
estimation will be more demanding. These parameters will also have some degeneracy with the parameters describing the background cosmological 
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dynamics. A joint Fisher analysis of MG parameters and cosmological parameters will require a dedicated study and will be presented elsewhere. 
As the Planck observations provide an accurate baseline standard cosmological model, using Planck prior may be a useful practical solution. 

The comparison of the linear and non-linear curves in Figures-fT4l and 1151 suggests that the non-linear contributions can make a significant 
effect. However, in the non-linear regime the covariance matrix should include the bispectrum and the trispectrum of the density field, which would 
reduce the Therefore, we can expect that a fully non-linear analysis, with a better modelling of the covariance matrix (e.g., from numerical 
simulations) would lower down the non-linear curves and make them closer to the linear result. However, because the departures from the ACDM 
power spectrum increase on weakly non-linear scales, non-linearities should still improve the discriminatory power of the analysis of 3D clustering 
as compared with the linear result. Here we follow a conservative approach as we estimate the constraints on the modified-gravity parameters by 
using the linear-theory 

We have studied the constraints in both f{R) and Dilaton models with and without redshift space distortion. We find that redshift space 
distortions only lead to a small broadening of the constraints. Below galaxy number density of N = 10“'*Mpc“^ the surveys rapidly lose their 
ability to discriminate. On the other hand, increasing the number density beyond N = 10“®Mpc“® does not lead to drastic improvements of the 
results. 

Our constraints are derived for a spectroscopic survey. Inclusion of photometric redshift error will degrade the discriminating power of the 
survey. 

Our constraints are based on lower order tangential modes I < 25. We find that the information contents in different t are highly degenerate. In 
the range 0.01 < k < 0.2h Mpc“^, we find that for more than three to four (logarithmic) bins in k the covariance matrix can become singular. 

Finally, our results are based on an uniform distribution of noise (constant N) while real surveys may have more complicated variations in the 
average galaxy number density which may depend on the radial and angular coordinates in the sky. 


6 DISCUSSION & FUTURE PROSPECTS 

In this paper we have studied the possibility of constraining modified gravity theories using galaxy clustering. We have studied two different 
modified gravity models: f{R) theories and Dilaton theories. In both cases we have used a specific parametrization and computed the as a 
criterion for constraining the departure from the ACDM models. We used the sFB transforms and the resulting 3D power spectrum Ce{k) for a 
range of £ and k values to constrain the model. We assume an all-sky coverage and a spectroscopic survey with a Gaussian selection function 
i/9(r) oc exp(—r^/ro), ro — 150/i“^Mpc. We fix the number density of galaxies to be N = 10“"^ Mpc“®. We find that the low ^ < 25 modes of 
C| {k, k') (with radial modes restricted to A: < 0.2 h Mpc“^) can constraint the parameter fug at a level of 2 x 10“® (3 x 10“®) with 3a confidence 
for n = 1(2). For the Dilaton models some of the parameters (/3o, mo) can be well constrained using galaxy clustering though there are others (s, r) 
which remain poorly constrained. The parametrization used by us depends on a tomographic approach. Combining constraints from higher ^ > 25 
modes can further reduce the error-bar and thus in principle make cosmological probes of gravity competitive with solar system tests. However this 
will require an accurate modelling of non-linear clustering in redshift space as well as the covariance matrix of the sFB power spectra. Our results 
are based on the linear power spectrum. However, we find that inclusion of nonlinear effect can drastically improve the x^- However we would like 
to emphasise that the contribution from higher order moments notably the trispectrum is not included in our covariance matrix. Although, this may 
not play an imprortant role in the quasi-linear regime, in the nonlinear regime such contributions can no longer be ignored. A detailed analysis will 
be presented elsewhere. Partial sky coverage, boundary conditions from specific survey geometry will also mean that in a realistic survey the sFB 
modes will have to be discrete and not continuous even in the radial direction. The average number density of galaxies and hence the noise too will 
be a function of the angular position on the surface of the sky. A pseudo-C^ (PCL) based approach can be useful in this direction. The errors due to 
photometric redshift determination can be readily included in our calculation. 

We have only considered two sets of MG theories. However, the results derived are quite general and can readily be extended to other MG 
scenarios such as the K-mouflage models or the Symmetron models. In addition, the sFB power spectrum can also be useful in constraining massive 
neutrinos, warm dark matter candidates or axionic dark matter through their footprints on the matter power spectrum as a function of redshift. 

The number of parameters required to specify a MG model can be high. Non-parametric techniques such as Principal Component Analysis 
(PCA) can be efficient to investigate linear combinations (ordered according to decreasing signal-to-noise) of parameters that can be extracted using 
a specific survey strategy. Extending the x^ based approach presented here, a more extensive analysis of survey optimisation covering a range of 
survey parameter and associated Fisher matrix analysis detailing parameter degeneracies will be presented elsewhere. In future we also plan to use 
a Bayesian model selection approach based on the evidences for acceptance or rejection of specific models. 

Physics of galaxy formation may be different in modified gravity models. Even in the ACDM cosmology, the bias associated with the galaxies 
may depend on the galaxy type and their environment, and may be stochastic. Such complications can only be dealt with numerical simulations. 
However, one may expect that at very large scale such effects may be less important. 

CMB mostly probes high redshifts where any modification from GR is expected to be sub-dominant. However, CMB is al so sensitive to modifi¬ 
cation s of gravity through secondary effects such as the Integrated Sachs-Wolfe (ISW) effect and the w eak lensing distortio ns jPlanck CollaborationI 
l2015al) . However the constraints are rather weak compared to what can be achieved by galaxy surveys jMunshi et al.ll2014h . 
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Weak lensing i s a very exciting possibility for constraining any departure from GR especially because it is free from the problem of galaxy bias 
iMunshi et alfcOOsh . However, in recent years many possible systematics have been discussed, e.g. from intrinsic ellipticity correlations or possible 
modifications of the power spectrum due to AGN feedback, gas cooling and stell ar feedback or neutrino physics (Mu nshi et al. 2015; in preparation). 
Such contaminations may induce appreciable bias in parameter estimation (e.g. lOsato, Shirasaki. & Yoshid3l2015h ). Despite all these systematics, 

weak lensing and galaxy clustering remain the most powerful probe of possible departures from G R on cosmological scales. _ 

Statistics of Lyman-a absorbtion have also been investigated in the context of f{R) gravity lAmold, Puchwein & SDrineelll2015h . Using cos¬ 
mological hydrodynamical simulations of f{R) gravity which include the flux probability distribution functions and the flux power-spectra, and an 
analysis of the column density and line-width distributions, as well as the matter power spectrum. It was found that Lyman-a statistics is rather 
insensitive to modification to gravity in f{R) models. Hence no competitive constraints are achievable using current data. Moreover, contamination 
by baryonic physics, associated with star formation and coolin g processes, imp ly that a very accurate modelling of such ingredients is required. 
Constraints from 21cm intensity mapping are more encouraging jHall et aD2013h . 

Observations by SDSS gal axy surveys in th e low to moderate redshift range 0.15 < 2 : < 0.67 is already being used to set interesting constraints 
on any departure from GR, e.g. lSel et alJ l l2014l) obtained fug < 4.6 x 10“® at the 95% confidence level. Euclid with survey characteristics similar 
to SDSS will probe deeper and wider parts of the sky thus improving the constraints by orders of magnitude. 
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